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ABSTRACT 


Prevention of pogo by means of an accumulator at the pump inlet is facilitated by a 
simplified analytical model that expresses the net in-phase force feedback in terms of 
the admittance of the accumulator, thereby providing a frame of reference for the use of 
judgment in the design. Treating the effects of the fuel and oxidizer systems as inde- 
pendent and additive and reckoning force feedback at the structure modal frequency only 
permit a simple description of the primary mechanism. A novel feature of the model is 
that the effect of outflow perturbations on tank -bottom pressure is included in the net- 
thrust function. 
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SUMMARY 

Prevention of pogo by means of an accumulator at the pump inlet is facilitated by a 
simplified analytical model that expresses the net in-phase force feedback in terms of 
the admittance of the accumulator, thereby providing a frame of reference for the use of 
judgment in the design. Treating the effects of the fuel and oxidizer systems as inde- 
pendent and additive and reckoning force feedback at the structure modal frequency only 
permit a simple description of the primary mechanism. A novel feature of the model is 
that the effect of outflow perturbations on tank -bottom pressure is included in the net- 
thrust function. 


INTRODUCTION 

The type of instability that has come to be known as "pogo” consists of longitudinal 
oscillation occurring near the end of the thrust period in liquid-propelled rocket vehicles. 
The mechanism is a feedback loop in which propellant -flow oscillations driven by the 
structure produce thrust oscillations that, in turn, drive the structure. The pogo prob- 
lem was studied in connection with the Gemini launch vehicle and satisfactory corrective 
devices were developed (refs. 1 and 2). Analytical models have been developed both for 
specific vehicles and for more general application, including a general algebraic model 
(ref. 3) and a general computer model (refs. 4 and 5). 

Reported pogo occurrences have the following characteristics in common: 

(1) The structure oscillates in the lowest -frequency mode in which the oscillation 
velocities of the engine and payload are opposite in phase, and the frequency of oscil- 
lation closely approximates the modal frequency. 

(2) Pogo occurs during the last part of the thrust period, usually shortly before 
cutoff. 



(3) Occurrence of pogo is associated with a condition of resonance or a near approach 
to resonance in the pump feed system of at least one propellant. 

Various methods of preventing pogo have been proposed, including 

(1) Increasing structural damping to make the system stable in spite of thrust feed- 

back 

(2) Changing ullage pressure to change pump gain and feed-system resonant fre- 

quency 

(3) Increasing feedline area to increase pressure force at the pump inlet 

(4) Use of a compliant device (accumulator) at the pump inlet to lower feedline re- 

sonant frequency 

(5) Use of a volume -compensating device to reduce or eliminate input from the 

structure to the propellant system 

(6) Use of compliant inserts or gas bubbles in the feedline to increase damping and 

lower resonant frequency 

Of these methods, the pump -inlet accumulator is probably the most practical and econ- 
omical. In its application the device has been designed to lower the fundamental re- 
sonance of the feed system to below the structure frequency and keep the second feed- 
system frequency above that of the structure throughout the thrust period. This method 
is satisfactory but may lead to overconservative requirements. 

The model presented herein is the outcome of a quest for an analytical representa- 
tion of pogo that would be simple enough to provide physical insight into the basic mech- 
anism and realistic enough to be useful as an aid in establishing design requirements 
for accumulators . Exploratory calculations using a bipropellant model similar to that 
of reference 3 indicated that the feedback effects of the two propellant systems are 
roughly independent and additive. It became apparent that the model could be simplified 
greatly by neglecting the interaction between propellant systems, treating each as if there 
were no oscillation in the other. A simple model based on this assumption was developed 
in the hope that such a model, requiring a minimum of calculation and providing a sim- 
ple description of the primary effect of accumulator parameters on propulsion feedback, 
would facilitate the selection of tentative design requirements and simplify the trial-and- 
error procedures required for optimization. 

The model contains two major simplifying assumptions in addition to those usually 
found in pogo models: that the effects of the two propellant systems are independent and 
additive, and that stability at any time can be determined by calculating force feedback 
only at the modal frequency of the structure. Otherwise, this model is very similar to 
that of reference 3 except that the accumulator admittance is included explicitly and the 
equations are cast in a form designed to reveal the effect of accumulator admittance on 
net force feedback. A new feature in the present model is the inclusion of the contribu- 
tion of the tank -bottom pressure oscillation to the retarding force exerted by the feed 
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system on the structure. 

The model is intended primarily for use as a design tool. Precise stability analysis 
requires the use of a complete bipropellant model including the propellant interaction 
effect. 


ANALYTICAL MODEL 
Assumptions 

The following assumptions are used in the development of the model: 

(1) The problem may be treated as linear for small perturbations about steady-state 
values . 

(2) Vehicle and propulsion -system parameters vary slowly with time and can be con- 
sidered to be constant for sinusoidal-perturbation analysis at any given time. 

(3) Oscillatory velocity and displacement of parts of the structure and flow pertur- 
bations in inlet lines are defined with respect to a reference system moving with the 
smoothed velocity and acceleration of the center of mass. Flow perturbations down- 
stream from the pump inlet are defined with respect to the conduit. 

(4) Compressibility is negligible between the pump inlet and the injector. 

(5) Friction losses in inlet lines are negligible. 

(6) Effective coupling between structure and propellant feed system occurs only 
upstream of the pump. 

(7) Pump-inlet cavitation behaves as a pure compliance. 

(8) In the pump equation, gain (m + 1) is real, and resistance R is independent of 
frequency. (Symbols are defined in appendix A.) 

(9) An adequate structure model exists, from which modal frequency co n , effective 
mass at gimbal point M n , and mode-shape factors cp n for selected structural elements 
may be calculated as a function of time for modes of interest. The structure model 
includes liquid in the propellant tanks, but the feedlines are empty. The lines are as- 
sumed to be closed at the tank bottom. Therefore, forces exerted by the liquid on the 
structure when the outflow oscillation is zero are included in the structure model. 

(10) Forces applied to the structure by the propellant system do not modify the mode 
shapes. 

In addition to the foregoing assumptions, which are common to most pogo models 
(refs. 3 and 4), the following special assumptions are used in the present model: 

(1) The component of the tank-bottom pressure oscillation caused by propellant out- 
flow oscillation has a negligible effect on the propellant-flow response to structure 
motion, but the resulting component of pressure force on the structure (not included in 
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the structure model) constitutes a significant part of the total reaction of the propellant 
system on the structure . 

(2) Effects of the two propellant systems are independent and additive. Oscillation in 
either may be determined by assuming no oscillation in the other, and the sum of forces 
on the structure caused by both independently is not substantially different from the 
effect of both acting together. 

This assumption permits a simple formulation of the basic feedback mechanism unob- 
scured by confusing second-order effects. It is justifiable because the determination of 
requirements for a corrective device requires that attention be focused on the primary 
effect of one propellant system. 

(3) The structure is massive and very lightly damped; therefore, the thrust- 
feedback function calculated at the modal frequency and only at this frequency may be 
used to determine stability. (This assumption represents a departure from the point of 
view of control theory because frequency is not regarded as an independent variable . ) 

The logical basis for this assumption is more fully developed in the following section. 


Stability Criterion 

A condition of neutral stability exists if, at any frequency, the mechanical impedance 
of the structure at the gimbal point equals the net -thrust transfer function. The mech- 
anical impedance is the complex ratio of driving force to velocity required for 
steady-state sinusoidal oscillation. It may be calculated as a function of frequency from 
the structure model. The net -thrust function Tj^/Xg is made up of the engine thrust 
oscillation minus the effect of forces applied to the structure as a result of flow oscil- 
lations in the propellant lines. It is expressed as the ratio of effective force to velocity 
at the gimbal point. 

The general nature of the variation in magnitude of these transfer functions with 
frequency at a particular time is illustrated in figure 1(a). The mechanical impedance 
is characterized by a very sharp minimum at the natural frequency of each structural 
mode. For bipropellant systems the net-thrust function generally has not more than 
two peaks, and these are broad compared to the narrow -band minima of Z . Since 
iT-^/XgJ must exceed j j for instability, it is clear that oscillation is possible only 
at frequencies very close to the natural frequency of a structural mode (modal fre- 
quency) . Moreover, the minimum value of Z m is ordinarily least at the lowest- 
frequency mode; therefore, the higher modes are less likely to be unstable. 

A comparison of the net -thrust and mechanical -impedance functions in the neigh- 
borhood of a modal frequency w n is shown in more detail in figure 1(b), in which the 
real and imaginary parts of both transfer functions are shown as a function of u>/a> n - 


4 



I 



Angular frequency, u) 


(a) General nature of the variation. 



Ratio of angular frequency to angular modal frequency, w/a> n 
(b) Detailed comparison. 

Figure 1. - Comparison of net -thrust transfer function and gimbal-point 
mechanical impedance. 


The criterion for stability is that the system is stable if the real part of T N /x g is less 
than the real part of Z m at the frequency at which the imaginary parts are equal. Be- 
cause of the very rapid change in the imaginary part of Z m with frequency (a conse- 
quence of large mass and very light damping) , only a small error is introduced if the 
comparison of the real parts of Z m and T N /x g is made at the modal frequency a> n 
instead of the frequency of phase coincidence. This is an important simplification be- 
cause it eliminates the use of frequency as a scanning variable. The net-thrust function 
at w n may be calculated as a function of a single independent variable, time. 

The mechanical impedance at the modal frequency is 
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The frequency and modal mass M n (with respect to the gimbal point) are obtained from 
the spring-mass model of the structure, but the damping factor £ n must be obtained 
from test data or estimated from experience with similar vehicles. Because of uncer- 
tainty in £ , it is convenient to present the results of stability analyses in terms of the 
damping factor required for neutral stability Zg. 

On the basis of the approximate stability criterion, ^ is given by 




C " 



2w M 
n n 


( 2 ) 


when T N /Xg is evaluated at the modal frequency o> n . The assumption that the effects 
of the two propellant systems are independent and additive permits the separation of the 
net -thrust function into two components 


T 


N 



(3) 


For either propellant acting alone, the thrust function is conveniently expressed in di- 
mensionless form in terms of two force -feedback parameters £ and v defined by the 
following equation: 


k + iv X (4) 

2w M x„ 
n n g 

In this equation F is the net effective force at the gimbal point resulting from oscil- 
lations in one propellant system at the frequency o> n . Since only the real part of the 
thrust function appears in equation (2) , the damping factor required for neutral stability 
is 


^f^o 


(5) 


Expressions for £ in terms of vehicle and propellant-system parameters are to be de- 
termined from an analysis of the propulsion system. 




Figure 2. - Schematic diagram of system. 


Inlet System 

A simple form of propellant system is illustrated in figure 2. It consists of a tank, 
a short feedline, a pump, a discharge line, and an injector. An accumulator is attached 
at the pump inlet for the purpose of stabilizing the system. The admittance of the accu- 
mulator is represented in general form by 


Y a = 


G + iB„ 


( 6 ) 


The thrust chamber is assumed to move with the gimbal-point oscillation velocity 
x (positive forward) . The structure oscillates at a modal frequency . The ratio of 
amplitude at any point to the gimbal-point amplitude is expressed by the mode-shape 
factor cp. In the model of the structure the liquid in the tank is represented by a single 
lumped mass having the oscillation velocity (p^g when the feedline is closed at the en- 
trance. With the feedline open, the flow perturbation entering the feedline is Q ; and, 

s 

therefore, the oscillation velocity of the liquid center of mass becomes 
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The pressure perturbation at the entrance to the feedline is, therefore, 


P t = ico n Ph cp q k g 


^Aq 

Vt : 


( 8 ) 


The second term is small compared to the first and may be neglected at this point, giving 


P t = ia, n ph< ? 



( 9 ) 


(It is usually assumed that the tank-bottom pressure perturbation is independent of the 
flow perturbation in the line (refs. 1 and 3) . Since Q g /S t « x g , this assumption is 
valid when the purpose is to determine the feed-system response to structure oscillation. 
However, if it is desired to determine the forces on the structure resulting from oscil- 
lation in the feed system, the flow -dependent term is important because the component of 
pump suction pressure in phase with produces a quadrature component in Q g that 
gives rise to a component of P^. in phase with and opposing the structure oscillation 
velocity.) 

Since the feedline is short compared to a/w n and is terminated in a low impedance 
at the tank, it may be treated as a lumped inertance given by 


L 


s 



( 10 ) 


The flow perturbation in the line (positive aft) is related to the pressure difference as 
follows 


Q = 

^s 


iw n pZ s 


< P t - P s> 


(ID 


The compliance of cavitation in the pump inlet is the rate of decrease in total 
volume of bubbles with increasing pump suction pressure. The response of volume to 



pressure is not purely compliant but is characterized also by an inertial effect arising 
from the acceleration of liquid about the bubbles and a lag resulting from heat transfer 
associated with evaporation and condensation. Since these effects cannot be estimated 
reliably and are thought to be minor, it is assumed that the cavitation behaves as a pure 
compliance, and the flow into the compliance is given by 

«b = *"n c b p s < 12 > 

An accumulator having an admittance Y (a complex function of frequency and op- 

a 

erating pressure) is attached at the pump inlet. The flow into the accumulator is 

«a - Y a P s < 13 > 

The flow perturbation entering the pump is defined with respect to the pump casing; 
therefore, the flow equivalent of the pump oscillation velocity is included in the equation 
for Qp. 


% ~ Q s + S sVg " ' y a (14) 

Substitution from equations (9), (11), (12), and (13) into equation (14) yields the low- 
pressure feed-system equation 


V 


= Ax, 


YP = S +-- tp )x 
s sl^p . ^q I g 



+ io> C, 
n b 



(15) 


This equation applies to the configuration of figure 2 when w Z /a < 0. 1 ir. In general, 

ns 

an analysis of any practical low-pressure feed system yields an expression of this form. 

For a longer line such as analyzed in reference 6 with 0. Itr < u> i /a < 0. 97 r, parameters 

n s 

A and Y are given by 



(16) 
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(17) 


Y - Y a + ‘"n C b + 




i 0)„pZ „ 
n s 


a 


For most systems A is real, and Y is imaginary for the unmodified feed system 
(damping negligible in inlet line) . With the accumulator attached, the combined admit- 
tance may be written as the sum of the accumulator admittance G & + iB & and the un- 
modified system admittance iB u> 


Y = G a + i(B a ♦ B u ) (18) 

Equation (15) shows that the system behaves as if the effect of structure motion were 
to inject into the feed system at the lower end of the line a flow oscillation equivalent to 
that produced by a piston of area A moving with the gimbal -point velocity x . The 
reaction of the feed system on the structure is, therefore, equivalent to a retarding force 
acting at the gimbal point given by 


F r = AP s (19) 

The expression AP includes the significant effect on the structure caused by the com- 
s 

ponent of tank-bottom pressure that results from propellant outflow oscillations. This 
fact is verified for the configuration of figure 2 as follows: Equations (8) and (11) are 
used to obtain an expression for P t in terms of x g and P g without using the approxi- 
mation of equation (9) . 


i“ ph Vg + “? l P S 


P t = 


S~ cp c 
s t n 




l S t cp t 


The above denominator is a scaler quantity that may vary only from 1. 0 to about 1. 1 in 
extreme cases; therefore, 


P t = iwph<p q x g 




I S 


t <Pt 
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The first term is included in the structure model, and the additional component of tank- 
bottom pressure due to propellant flow oscillations is given by 


h fs fa p 

' s t i’i 


and the effective propellant-system reaction at the gimbal point is 


F d = S 
R s s^p 


+ s t p^ t - S g 




<Pc 


= AP g 


In order to obtain expressions for the pump-flow and suction -pressure transfer 
functions before proceeding with the analysis of the discharge system, it is assumed at 
this point that the relation between pump throughflow and suction pressure may be rep- 
resented by an effective impedance (the assumption of independence of fuel and oxidizer 
systems is implied), 



( 20 ) 


Equations (15) and (20) describe a dividing network for which the ratio of output to input 
is given by 


% = 1 

Ax 1 + YZ’ 
& 

The pump-flow transfer function is 

Qp _ A 
x g 1 + YZ’ 

And the suction-pressure transfer function is 

^ _ AZ* 
x g 1 + YZ’ 


( 21 ) 


( 22 ) 


(23) 
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Pump and Engine 


The dynamic behavior of the pump is represented by an equation of the form (ref. 1) . 

P d = (“ + 1)P, - Z p Qp (24) 

Limited unpublished results from tests conducted at Lewis indicate that the real part of 
Zp may be frequency -dependent and that m + 1 may be complex. However, for the 
purpose of the present analysis, it is assumed that m + 1 is real and that Z p can be 
represented by 


z p = R p + iuL p (25) 

In the absence of dynamic test data, m + 1 and Rp may be estimated from steady -state 
performance data (see ref. 1), and a rough approximation of L p may be calculated from 
the internal geometry of the pump . 

Since the injector resistance is small compared to the characteristic impedance of 
the discharge line, it is assumed that the system downstream of the pump is incompress- 
ible; therefore, Qp represents also the flow perturbation at the injector. The discharge 
line and injector are represented by a lumped impedance in the pressure equation: 

P d = P c + Z d% < 26 ) 


where 


Z d - R d + iwL d 

Some uncertainty exists concerning line resistance to oscillatory flow when the mean 
throughflow is in the turbulent regime. Linearization of the steady-state pressure-flow 
relation is valid only if the period of oscillation is long compared with the time required 
to establish the steady-state turbulent -flow velocity profile. 

An equation for the pump-inlet pressure perturbation is obtained from equations (24) 
and (26) 


Z + Z 

P =_£ do 

s m + 1 ^ 


m + 1 


(27) 
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The pressure oscillation in the thrust chamber results from flow oscillations of fuel and 
oxidizer and is characterized by a time lag that is small compared to l/w n - The average 
combustion delay time is added to the first-order lag of the chamber volume to obtain 
the time constant r . The pressure oscillation is given by 


p _ H fipf + gogpo 


(28) 


The partial-resistance coefficients H Q and H f are defined as the zero -frequency 
variation of chamber pressure with propellant flow for one propellant, assuming steady 
flow of the other. These coefficients include the effect of variation in oxidant-fuel ratio 
on characteristic velocity. They may be determined from engine performance data as 
described in reference 3. (In the notation of ref. 3, H is represented by C*p/A^g.) 

In order to calculate an effective impedance for each propellant looking from the in- 
jector into the thrust chamber, the interaction implied in equation (28) is neglected, and 

2 2 

oscillation in only one propellant system is assumed. Noting also that w n r c « 1, the 
effective impedance is, for either propellant, 

— = H(1 - iw r ) (29) 

Qp 

The effective pump-inlet impedance z' required for use in equations (22) and. (23) is 
obtained from equations (27) and (29) 

— = Z* = R’ + iX* (30) 

% 

where 

R’ = -A- (R + R d + H) 
m + 1 p a 


X f =-^(L + L d - H r c ) 
m + 1 • 


13 



Net-Thrust Transfer Function 


The transfer function relating chamber pressure to gimbal -point oscillation velocity 
is obtained by assuming that the effects of the two propellants are additive, though each 
is determined as if the other were zero. Using equation (22) in equation (28), 


^c = 

HA 

+ 

f 

HA 1 

• 

g 

(1 + iwr c )(l + YZ’). 

.(1 + ia>r c )(l + YZ')] 

l 


(31) 


The driving force on the structure at the gimbal point produced by the chamber -pressure 
oscillation is S E P c’ where is the ratio of mean thrust to mean chamber pressure. 
Because the reaction of each feed system on the structure is equivalent to a retarding 
force AP g at the gimbal point (eq. (19)), the net -thrust perturbation is 


T N " S E P c 


A f P f - A P 
f sf o so 


(32) 


The chamber -pressure and suction-pressure transfer functions (eqs. (31) and (23)) 
are used in equation (32) to obtain the net -thrust transfer function 



s e ah a 2 z' 1 

+ 

f 

s e ah a 2 z’ 1 

k e 

.(! + i co n T c )( 1 + YZ’) 1 + YZ'_ 

- (1 + ico n T c )(1 + YZl) i+YZ’J 


(33) 


This equation is in the same form as equation (3); each term represents the thrust 
function for one propellant. In order to show explicitly the effect of an accumulator, 
equation (18) is used for Y in equation (33), and the dimensionless force-feedback 
function for one propellant (eq. (4)) is written in a form intended to reveal the effect of 
G and B . For either propellant, 

cl a 


£ + iv = 


V H f 1 ' 1 1 

2w n M n V 1 + lw n T c 



R’ 




, x * - ( B a + b u ) i z *r 

l-l 


R’ + G |z'| 2 

cl 


(34) 


14 



In order to facilitate the intelligent selection of G and B , equation (34) is trans- 

a a 7 

formed in a manner intended to reveal more clearly the effect of various system param- 
eters 


i+i^fie^ + ^cos* 
KX 


(35) 


This form was obtained by representing the four factors in equation (34) as follows: 


s e ah 

2w n M n 


= K 


(36) 


R 


R* + gIz ’| 2 


= X 


‘-‘I 


A|z 'i =ne ie 


1 + ‘"n T c s e hr ' 


(37) 


(38) 


1 


1 - 



-(B a + B u>l 2 ’| 2 

R' + G a |z '| 2 


1 Up . 

= e ^ cos xp 

1 - i tan xp 


(39) 


The parameter K is a scale factor that expresses, in units of £, the gross thrust for 
direct in-phase coupling as would be realized with undamped inlet -system resonance, 
pure resistance on the discharge side, and no lag. The factor A describes the reduc- 
tion in feedback amplitude caused by damping in the accumulator. The complex factor 
i 0 

fie represents, in polar form, the combined effects of the phase shift caused by 
thrust-chamber lag and discharge-system reactance and the relative back force exerted 
by the inlet system. The effect of inlet-system resonance is described in terms of the 
detuning angle \p. 

Representation of the normalized force-feedback function (| + iv)/KX in the com- 
plex plane (fig. 3) illustrates the effect of fi, 9, and \f/. Parameters fi and 6 are 
independent of resonance conditions in the feed system. Therefore, a single point rep- 

; n 

resents fie at a given time for all possible combinations of B„ and G . The locus 

a cl 
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Figure 3. - Complex-plane representation of normalized force- 
fee&adc function (£+ iv)/KA. 

of (| + iy)/KX is a circle constructed on the radius vector £2e iy as a diameter. For 
large negative values of B & + B u , the angle i p approaches 90°. With increasing 
B a + B u , \{/ decreases, passing through zero when B a + B u - x ’/| Z f |^ and approaching 
-90° at large positive values of B & + B u . The horizontal extent of the circle defines the 
limits within which |/KX may vary with changes in inlet-system admittance, and since 
X cannot exceed 1, the limits of £ are 

«max = 8 + 1) < 40 > 

Wn'f^ 089 - 1 ) < 41 > 

The variation of £ between these limits is obtained by taking the real part of equa- 
tion (35) . 

ij = KXO cos (6 + xjJ) cos (42) 
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APPLICATION 


Recall that the structure damping factor required for neutral stability is obtained by 
evaluating | for both propellants at the structure modal frequency o> n . 

C C - «£ * «o < 5 > 

Since the accumulator conductance G„ and susceptance are instrumental in deter- 

a a 

mining if/ and X, the variation of | between the limits given by equations (40) and (41) 
can be influenced by specifying G & and B & , and £ may be confined within allowable 
limits. 


Effect of Accumulator on Thrust-Feedback Function 

The variation of £ with the combined inlet-system admittance Y = G. + i(B„ + B„) 

can be understood by means of the complex -plane diagram (fig. 3) . For the unmodified 

system, Y = B u and X = 1. As B u increases from large negative values, the detuning 

angle if/ decreases from near 90° and the point representing (£ + iy)KA moves clockwise 

along the circle. The maximum value of £/KX occurs when \f/ = -9/2, and the minimum 

when ip = -90° - (0/2). The variation of £ between its limits depends on the variation 

with time of the modal frequency c o„ in relation to the inlet -system frequency u>„ (at 

which B = 0). Two examples are shown in figure 4. Figure 4(a) represents a short 

feedline with pump-inlet cavitation increasing near the end of firing. The inlet-system 

frequency decreases, crossing the structure frequency. In this case £ starts small and 

positive, increases to the upper limit as the frequencies cross, decreases quickly to the 

lower limit, and rises slowly toward zero. Figure 4(b) represents a long feedline with 

decreasing cavitation (increasing P_) . The sequence of events is reversed as oj 

s s 

crosses a> n from below. 

Also shown in figure 4 is the effect of small undamped and lightly damped accumula- 
tors. In both cases the resonant frequency of the accumulator is higher than o> n , making 
B positive at all times. In figure 4(a), the maximum of £ occurs earlier in flight as 

a 

a result of lowering the inlet-system frequency a>_. In figure 4(b), the pattern of vari- 

s 

ation in £ is shifted in the opposite direction. With damping, the peak value of £ is 
below the limit given by equation (40) , and the system passes through resonance more 
slowly. 

If the inlet-system frequency tends to track the structure frequency as in figure 4(b), 
a small undamped accumulator is sufficient to provide adequate frequency separation at 
all times. The system represented by figure 4(a) poses a more difficult problem because 
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Figure 4. - Effect of accumulator on thrust feedback. 


a very large accumulator would be required to provide frequency separation at lift-off. 

The situation is not as bad as it looks, however, because is much higher than co„ 

only if the feedline is short. And a short line implies a large value of A in the early 

part of the thrust period. In these circumstances the back force AP_ is large, and 

| „„ is so small that resonance may be allowed during the early part of the thrust 
max 

period. In this case damping in the accumulator is advantageous because it decreases 
the peak value of £ • 


Accumulator Design Considerations 

Because of the wide variety of possible relations among several parameters, it is 
not feasible to propose a general procedure for selecting design conditions. Optimal 
design requires judgment, experience, and at least a minimum of trial -and-error pro- 
cedure. The present model is intended to aid in the design procedure in two ways: (1) by 
providing a framework within which the designer may exercise judgment in selecting a 
design point, and (2) by providing a means of verifying tentative designs by means of 
relatively simple calculations. 

The pattern of in-phase force feedback for one propellant system is defined by the 
curve of £ as a function of time . Addition of an accumulator to the system gives rise 
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to a different curve, but all such curves are included between the envelopes defined by 

£ and £ . . Examination of the envelope curves and the unmodified £ curve for 

^max ^min 

both propellant systems provides the information required to decide whether either or 
both systems require modification. 

The curves of £ in figure 4 show the general nature of changes caused by increasing 
B. These curves and a study of the relation between figures 3 and 4 should suggest the 
character of changes in £ to be expected for systems different from those illustrated 
in figure 4. Based on a qualitative understanding of the shape of realizable £ curves, 
a single critical point can be selected for design. Examples of the type of design con- 
dition that may be specified are the latest time at which £ is allowed to reach the upper 
limit or the time at which £ is zero. Either of these conditions defines a limiting 
value of B„ + B„ at the selected time . The design condition also may be specified as a 

a U. 

limit less than £ for a peak value of £ at a given time, defining also a minimum 
^ max 

required value of G„ . 

Based on the design condition, the accumulator configuration and tentative values 
of the circuit elements are selected. Various configurations are possible, and a variety 
of admittance -frequency characteristics are realizable by combining compliance C, 
inertance L, and resistance R in various ways. The accumulator conductance G a 
and susceptance B are calculated for the equivalent circuit by methods applicable to 
electrical networks. Some possible configurations and the relation between the ideal 
elements R, L, and C and their physical counterparts are discussed in appendix B. 

The tentative design is used to calculate G„ and B as a function of time, taking 
into account variations of circuit elements and frequency w n with time, and these 
values are used to generate a curve of £ for the modified system. Additional checkout 
calculations can be made by using upper and lower limits for uncertain parameters. 

Similar procedures may be applied to the second or higher structural modes if 
doubt exists as to their stability. The fact that higher-mode oscillation has not been 
observed suggests the possibility that £ is generally less than i? except in the 
lowest mode. 

A specific example, illustrating some possible design procedures, is presented in 
appendix C. 


CONCLUDING REMARKS 

Interaction between fuel and oxidizer systems (due to the thrust chamber being 
common to both) has been neglected in the present model, as has the frequency shift 
caused by the quadrature component of net-thrust feedback. By means of these approx- 
imations, the model is reduced to a simple form that reveals the basic feedback mechan- 
ism unobscured by second-order effects, and indicates clearly the changes in thrust 
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feedback that may be realized by means of pump-inlet accumulators. This model includes 
also an improved method of reckoning the effective force exerted on the structure by the 
propellant systems, eliminating the error resulting from the assumption (in previous 
models) that tank-bottom pressure is independent of perturbations in outlet flow. 

Because of the simplifying assumptions used in its development and the difficulty 
anticipated in determining critical parameters such as pump gain and compliance, the 
model is not suitable for predicting instability in cases that may be marginal. However, 
this model was developed with a view to a practical solution of the longitudinal-oscillation 
problem. From this viewpoint an exact determination of stability limits is not required 
because prevention of instability (pogo), rather than prediction, is the aim. If prevention 
of pogo is deemed of sufficient importance to warrant modification of the propellant 
systems, the model may be used to define the required admittance of pump-inlet accu- 
mulators for fuel or oxidizer or both. The model may be used also to check the stability 
of the modified system, since prudent design requires that the in-phase thrust feedback 
be so small that moderate percentage errors in its determination would be tolerable. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 2, 1969, 

120-27-04-27-22. 
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APPENDIX A 


SYMBOLS 


A 

2 2 

input coefficient, in. ; m 

R’ 

effective resistance looking into pump 

a 

sonic speed, in. /sec; m/sec 

r 

radius, in. ; m 

B 

susceptance, in. ^/lb-sec; 

S 

2 2 
area, in. ; m 


m VN-sec 

t n 

net -thrust perturbation, lb; N 

C 

compliance, in.Vlb; m^/N 

V 

3 3 

volume, in. ; m 

F 

force, lb; N 

V 

velocity, in. /sec; m/sec 

G 

conductance, in. ^/Lb-sec; 
m 5 /N-sec 

X 

x’ 

5 5 

reactance, lb-sec/in. ; N-sec/m 
effective reactance looking into pump, 

H 

thrust -chamber partial - 
resistance coefficient, 


lb-sec/in. N-sec/m^ 


c r 

lb-sec/in. ; N-sec/m 

X 

displacement, in.; m 

h 

liquid level, in. ; m 

X 

velocity, in. /sec; m/sec 

i 

V 7 ! 

X 

acceleration, in./sec^; m/sec^ 

K 

force-feedback scale factor, 

Y 

admittance, in./lb-sec; m°/N-sec 


nondimensional 

Z 

impedance, lb-sec/in. N-sec/m^ 

k 

kinematic factor (see ref. 6), 
nondimensional 

z’ 

effective impedance looking into pump, 
lb-sec/in. *; N-sec/m^ 

L 

ft r A 

inertance, lb-sec /in. ; kg/m’ 

y 

ratio of specific heats 

l 

length, in. ; m 


damping factor, nondimensional 

l' 

effective length, in. ; m 

0 

force -feedback phase angle at re- 

M 

effective (modal) mass, 


sonance, deg; rad 


O 

lb-sec /in.; kg 

K 

spring rate, lb/in. ; N/m 

m + 1 

P 

pump gain, nondimensional 
2 

pressure, psi; N/m 

X 

force -feedback amplitude reduction 
factor, non dir" isional 

Q 

q 

volume flow rate, in. /sec; 
m Vsec 

V 

quadrature component of force- 
feedback function, nondimensional 

R 

resistance, lb-sec/in. ®; 
N-sec/m^ 


in-phase component of force -feedback 
function, nondimensional 
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p 

2 4 3 

density, lb -sec /in ; kg/m 

h 

thermal 

a 

pressure-loss coefficient, 

9m 

imaginary part 


nondimensional 

m 

mechanical 

T 

time constant, sec 

n 

modal 

<P 

mode -shape factor (gimbal 
point = 1.0), nondimensional 

nl 

nonlinear 


o 

oxidizer 

* 

force -feedback detuning angle, 




deg; rad 

P 

pump 

to 

force -feedback amplitude factor, 

q 

liquid center of mass 


nondimensional 

R 

propellant system reaction 

w 

angular frequency, sec”* 

Re 

real part 

Subscripts: 

s 

suction line or pump inlet station 

a 

accumulator 

t 

tank 

b 

cavitation 

u 

unmodified 

C 

critical (stability boundary) 

w 

working 

c 

thrust chamber 

1,2,3 

circuit elements 

d 

discharge system 



f 

fuel system 



S 

gimbal point 
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APPENDIX B 


SOME PRACTICAL ASPECTS OF ACCUMULATOR DESIGN 
Accumulator Configurations 

Perhaps the simplest form of accumulator consists of a volume of noncondensable 
gas connected to the propellant line by a short liquid-filled duct, as shown in figure 5. 

If the propellant is a liquid storable at ambient temperature and in which air is not highly 
soluble, the gas charge may consist of the air entrapped when the system is filled 
(fig. 5(a)). Otherwise, gas may be drawn from existing storage on the vehicle or from a 
local storage tank, as shown in figure 5(b) . Storage within the accumulator in a flexible 
container, as in figure 5(c), is another possibility. If the gas is initially at atmospheric 
pressure, the entrapped volume must be large enough to allow for compression to oper- 
ating pressure. Moreover, the possibility of some loss of gas during the starting trans- 
ient must be examined. 



(a) Entrapped air. 



Be I lows 7 


Gas ■ 


(b) Gas storage tank. (c) Flexible container. 


Figure 5. - Accumulator configurations. 

Relatively larger working volume may be obtained by using the arrangement illus- 
trated in figure 5(b), in which the gas charge is stored under pressure in a separate 
chamber. The connecting valve is opened after operating pressure is established. The 
resistance and inertial reactance of the connecting tube and valve can readily be made 
negligible compared to the compliant reactance of the gas remaining in the storage cham- 
ber. Thereby the effective compliance is made equal to that of the total volume of gas 
in both chambers. If damping is desired, use of a porous plug in the connecting tube 
provides a convenient means of realizing approximately linear resistance. Use of the 
connecting tube as an inertance element is not considered to be practical because of the 
low density of the gas. If the pressure decreases with time, it may be desirable to con- 
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trol the liquid level at a point at least one radius above the top of the duct to allow for 
finite oscillations without getting gas in the main duct. Liquid-level control may be 
accomplished by the use of a small secondary duct large enough to permit the escape of 
excess gas but so small that its impedance is much greater than that of the main duct 
(fig. 5(b)). 

The foregoing are only examples of some devices that might be employed. Actual 
design requires familiarity with the particular propellant system and associated hard- 
ware. 


Physical Elements 

The parameters that must be realized within specified limits are the susceptance 

B and conductance G . These parameters can be determined as a function of time at 
a a 

operating pressure and structure modal frequency if the device can ne represented by a 
fluid circuit consisting of approximately linear inertance, compliance, and resistance. 
Various physical elements that may be employed are represented by corresponding com- 
binations of ideal circuit elements, as shown in figure 6. 


Ic 

(a) Adiabatic gas volume. 
L r 


ICj X c 2 


(b) Packed gas volume, 
o — nnnnp^ 


(c) Bellows or spring-loaded piston. 


(d) Liquid-filled duct. 


(e) Porous slab or plug. 


Figure 6. - Circuit representation of elements. 


Adiabatic gas volu me. - If the thermal time constant of the gas in the enclosure is 
much longer than the period of oscillation, the gas may be represented as a simple com- 
pliance proportional to the working volume 


C = — (Bl) 

Y P 

At finite oscillation levels the nonlinearity of the pressure -volume relation introduces 
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harmonic distortion, but the effective value of compliance at the fundamental frequency 
is not strongly dependent on amplitude. Compliance from gas compressibility always 
acts in the circuit as a shunt to ground. 

Packe d gas volume . - If special means are provided to promote heat transfer to 
and from the gas, as by packing the volume with a material such as steel wool or fine 
wire and making the heat capacity of the solid large compared to that of the gas, the dy- 
namic behavior is represented by the circuit of figure 6(b). The adiabatic compliance is 
represented by Cj, + Cg is the isothermal compliance, and RC 2 = is the thermal 
time constant. Damping is significant at the modal frequency if o> n is in the neighbor- 
hood of 1. This means of providing resistance is desirable because the admittance is 
approximately independent of the flow -oscillation amplitude. 

Bellow s or piston . - A bellows or spring-loaded piston is represented by a com - 
pliance and an inertance in series, as shown in figure 6(c). The compliance is 



where S is the effective area dV/dx and k the spring rate dF/dx. The inertance is 


where M is the effective mass (e. g. , the free end plus one -third of the convolutions of a 
bellows) . 

Liquid-fil led duct . - The inertance of the liquid between the compliance and the pro- 
pellant stream is probably the most difficult element to estimate accurately. For a 
straight circular duct of length l , cross-section area S, and radius r opening at both 
ends into large volumes of placid liquid, the inertance is 

L = (B4) 

S 

where V -l + 2A l and Al =0.85r. 

The end corrections A l represent the inertial effect of the liquid beyond the ends of 
the duct expressed as an equivalent extra length. The same equation applies approxi- 
mately to curved ducts, l being the centerline length, when the radius of curvature of 
the centerline is greater than about 2r. For tapered ducts, pi / S is replaced by 
/p/S d l. 

Where the duct joins the propellant system, the main stream moves past the junction 
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at high speed, and the influence of the strong sh<\„r layer on the end correction cannot be 
determined analytically. In the absence of experimental data for liquid systems, uncer- 
tainty is unavoidable, but some basis for a qualitative estimate is found in the results of 
a study of the frequency shift of Helmholtz-resonator-type sound absorbers with increas- 
ing flow velocity in an air duct (ref. 7). It was found that the effective inertance of holes 
in a thin plate decreases at high tangential velocities to about 37 percent of its theoretical 
value (1.7 r). If this result is interpreted as a change in only the outside AZ , it would 
indicate a value less than zero. Intuitively, it appears unreasonable, with a duct of 
finite length, to have a negative end correction. But the result suggests strongly that the 
high velocity past the opening greatly reduces the end correction. It is suggested, there- 
fore, that an end correction of zero be used for design calculations with a generous 
allowance for uncertainty. 

No definite rules can be given for estimating the inertance within the accumulator 
from the duct to the gas volume. If a large area change occurs, the standard end cor- 
rection (0.85 r) is a good approximation. Otherwise, an estimate of effective cross- 
sectional area as a function of axial length may be used as a basis for calculation. 

Resistance elements. - Since the flow in and out of the accumulator is entirely os- 
cillatory, the action of a resistance element having a square-law characteristic is 
strongly nonlinear, the effective resistance increasing with the first power of the flow 
amplitude. Because of the difficulty in predicting characteristics under operating con- 
ditions with variable resistance, and because of the frightening possibility that high 
effective resistance might in some circumstances prevent recovery from transient struc- 
tural oscillations, it is important that resistance elements be as nearly linear as prac- 
ticable . 

Linearity of resistance may be assured by using resistance elements in which lam - 
inar flow prevails at the largest flow -oscillation amplitude likely to result from extra- 
neous disturbances. A practical form of resistance element consists of a slab, plug, 
or sheet of porous material with a great number of small flow channels acting in parallel. 
Determination of the required combination of area, thickness, pore size, and flow-area 
fraction must be based on test data or reliable specifications for the material. Use of a 
porous -metal element as a linear resistance is discussed in reference 8. 


Duct Area Requirement 

Another problem that requires examination is the effect of nonlinear resistance 
caused by duct entrance and exit losses at high oscillation amplitude. If the duct contains 
a linear resistance element, the duct area will, in all probability, be large enough to make 
the nonlinear resistance negligible at any reasonable oscillation amplitude. However, if 
the accumulator is designed for undamped operation, or if a compound accumulator with 
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resistance in the gas-filled connecting tube or some other arrangement using an un- 
obstructed main duct is contemplated, the nonlinear resistance imposes a minimum 
area requirement. 

A relation between the effective nonlinear duct resistance R and the accumulator 
flow amplitude is obtained by applying the square -law pressure -drop relation to a sinus- 
oidal flow oscillation (see ref. 9). 


(b5> 

m 3ir g 2 

Experimental values of the pressure -loss coefficient a (reciprocal of the square of the 
discharge coefficient) for oscillating flow of water are given in reference 9. At peak 
oscillation velocities of over 0. 5 meter per second, the results show reasonably good 
agreement with accepted values for continuous flow. In particular, o is approximately 
2. 6 for an orifice in a thin plate and 1. 5 for length-diameter ratios greater than 1. 

The minimum required duct area depends on the maximum expected flow amplitude 
and the maximum allowable increase in resistance. 

A suggested criterion for allowable resistance is 


when 


Q = 



CO 


n 


where the sinusoidal acceleration x^ is the maximum expected quasi-sinusoidal re- 
sponse of the structure to extraneous disturbance (ringing level). 

Using these expressions in equation (B5) gives the minimum required area 


S = 



1.7 crpAx 
Z ’-n 


g 


(B6) 
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APPENDIX C 


EXAMPLE 

The general configuration of the vehicle selected as an example is shown in figure 7 . 
Hydrocarbon fuel is stored in the upper tank and liquid oxygen in the lower. Values of 
pertinent vehicle and propellant -system parameters are presented in table I. All di- 
mensional parameters are given in inch - pound force - second units. Equivalent units 
in the International System are listed in appendix A. (Since this is an hypothetical ex- 
ample for illustrative purposes only, the actual numbers are of no special significance. 
Therefore, use of both systems of units in text and and tables would serve no useful pur- 
pose.) 

The solution of any design problem is subject to certain constraints imposed by the 
practical aspects of the situation. It is assumed in the present case that the following 
conditions have been established: 
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TABLE I. - HYPOTHETICAL VEHICLE PARAMETERS 


(a) Fixed parameters 


System 

Density, 

Thrust- 

Suction- 

Suction- 

Sonic 

Total 

Total 

Modal 

Kinematic 

Ratio of 

Thrust- 



chamber 

line 

line 

speed, 

resistance, 

inertance, 

velocity 

factor, 

mean 

chamber 


lb-sec /in. 

partial- 

area, 

length, 

a s> 

R p + R d> 

L P + L d- 

ratio of 

k 

thrust to 

time 



resistance 

coefficient, 

S s’ 

in. 2 

l s’ 

in. 2 

in. /sec 

lb-sec/in. 

lb-sec^/in. ^ 

pump, 


mean 

chamber 

constant, 

V 



H, 








pressure, 

sec 



lb-sec/in. 








S E’ 

in. 2 


Oxygen 

1.06X10' 4 

0.034 

36 

30 


0.110 

9X10' 4 

1 

1.0 

0 

285 

0.0025 



Fuel 

.795 

.030 

22 

320 

2.8X10 4 

.165 

14.1 

1.0 

-.78 

285 

.0025 


(b) Time-dependent parameters 


Burn time, 
percent 

Modal 

angular 

frequency, 

n’ 

sec“* 

Effective 

modal 

mass, 

M , 
n’ 

9 

lb-sec/in. 

Modal velocity 
of liquid 
center of 
mass, 

Pump gain, 
m + 1 

Liquid level, 
h, 

in. 

Suction 

pressure, 

P s’ 

psi 

Caviation 

compliance, 

c b> 

in. 5 /lb 


Oxygen 

Fuel 

Oxygen 

Fuel 

Oxygen 

Fuel 

Oxygen 

Fuel 

Oxygen 

Fuel 



0 

72 

2600 

1 . 

0 

-0.50 

1.0 

1 . 

0 

250 

200 

61 

62 

0.04 

0.02 

60 

92 

500 



-.24 

1.3 



110 

83 

45 

42 

.20 

.04 

80 

104 

180 



.18 

1.8 



63 

50 

-43.5 

50 

.40 

.03 

95 

126 

50 



.50 

3.0 



27 

23 

41.5 

61 

.75 

.021 

100 

148 

25 


r 

.58 

5.4 


r 

15 

15 

40.5 

65 

1.50 

.019 


N> 

CO 



TABLE n. - PROPELLANT SYSTEM CALCULATIONS 


(a) Calculations based on equations (40) and (41) 


System 

Burn 

Area, 

Effective 

Effective 

Force- 

Driving-force 

Retarding- 

Force- 

Force- 

In-phase 


time, 

mT 2 

resistance 

reactance 

feedback 

term, 

force 

feedback 

feedback 

component 


percent 


looking 

looking 

scale 

i-i*l 

term, 

amplitude 

phase 

of force- 




into pump, 

into pump, 

factor, 

R‘ 

A 1 Z* I 2 

factor, 

angle at 

feedback 




R 



K 

1 + iwT 

s e hr 

a 

resonance, 

function, 




lb-sec/in. 5 

lb-sec/in. ® 


c 


6, 












deg 

Maximum 

Minimum 

Oxygen 

0 

336 

0.144 

0.059 

0.0087 

0.90 - i0.57 

5.85 

4.98 

-173.4 

0.00015 

-0.0431 


60 

168 

.111 

.058 

.0163 

0.84 - i0.72 

3.80 

3.04 

-166.5 

. 00067 

-.0488 


80 

112 

.080 

.047 

.0290 

0.79 - i0.79 

1.24 

.91 

-119.4 

.00675 

-.0198 


95 

68 

.048 

.034 

.0526 

0.71 - 10.94 

.56 

.96 

-81.2 

.0290 

-.0212 


100 

54 

.027 

1 .022 

.0705 

0.61 - il . 06 

.26 

1.12 

-71.5 

.0520 

-.0270 

Fuel 

0 

40 

0. 195 

0.096 

0.0009 

0.88 - i0.65 

1.13 

0.70 

-110.9 

0.00021 

-0.00043 


60 

42 



.123 

.0036 

0.81 - iO. 81 

1.33 

.96 

-122.5 

.00080 

-.00265 


80 

37 



.139 

.0084 

0.76 - iO. 91 

1.27 

1.04 

-119.0 

.00225 

-.0065 


95 

32 



.168 

.0215 

0.66 - il. 07 

1.26 

1.22 

-119.3 

.00671 

-.0196 


100 

31 



.197 

.0353 

0.55 - il. 00 

1.40 

1.31 

-130.3 

.00820 

-.0381 


(b) Calculations based on equation (42) 


System 

Burn 

Unmodified 

Force- 

Force- 

Force- 

In-phase 


time, 

susceptance, 

feedback 

feedback 

feedback 

component 


percent 

in. 5 /lb-sec 

detuning 
angle 
tangent, 
tan \p 

detuning 

angle, 

deg 

phase 
angle, 
\f/+ 9, 
deg 

of force - 
feedback 
function, 
$ 

Oxygen 

0 

-154 

26.4 

87.8 

-85.6 

0.00012 


60 

-105 

15.3 

86.3 

-80.2 

.00055 


80 

-67 

7.78 

82.7 

-36.7 

.00273 


95 

4.5 

.39 

21.2 

-60.0 

.0234 


100 

45 

-1.23 

-50.9 

-122.4 

-.0267 

Fuel 

0 

-7.7 

2.37 

67.1 

-43.8 

0.00018 


60 

-2.0 

1.18 

49.9 

-72.6 

. 00067 


80 

-.86 

.96 

44.0 

-75.0 

.00163 


95 

1.6 

.328 

18.2 

-101.1 

-.00480 


100 

4.0 

-.58 

-30.1 

-160.4 

-.0377 
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(1) The accumulator is to have the configuration shown in figure 5(c), using a flexible 
bellows precharged with nitrogen. 

(2) Friction in the duct and bellows is negligible, but damping may be provided, by 
packing the gas volume with a mesh of fine wire to promote heat transfer (see appendix B) . 

(3) The structure damping factor is estimated to be 0.015, and a safety factor of two 

is regarded as necessary; therefore, the maximum allowable value of is 0.0075. 

Envelopes of maximum and minimum values of the dimensionless in-phase thrust- 
feedback parameter £ were determined for both propellant systems by using equa- 
tions (40) and (41). These calculations are summarized in table n(a). Additional calcu- 
lations using equation (42) to determine £ for the unmodified systems are listed in 
table H(b) . Results are plotted as a function of percent burn time in figure 8. 

The curve of 4 0 for the unmodified oxygen system in figure 8(a) shows that the 
oxygen system alone would make the vehicle unstable at the 1. 5 percent level from 90 to 
96. 5 percent burn time. The corresponding curve for the fuel system (fig. 8(b)) shows a 
maximum of 0. 002 to 0. 0025 at about 85 percent burn time. These results indicate 
that the oxygen system requires modification; but the fuel system can be let alone since 
it contributes little to the instability. Therefore, the accumulator for the oxygen system 
is to be designed to stabilize the vehicle with the fuel system unmodified. In order to 
provide for consideration of the fuel system in selecting a design point for the oxygen 
system, a curve showing + £ was plotted on the oxygen-system diagram 

(fig. 8(a)) . This curve shows, for a portion of the flight, the combined effect of the un- 
modified fuel system plus the maximum possible effect of the oxygen system. 

Within the assumed "ground rules", a reasonable procedure would be to design an 
undamped accumulator to make the maximum system feedback equal to 0. 0075 and then 
calculate the effect of adding damping (and increased compliance) by heat conduction, in 
the hope of improving the margin. 

Circuits representing the specified accumulator types were obtained by combining 
the elements described in appendix B. The circuit of the undamped accumulator with 
bellows is shown in figure 9(a) and a simplified equivalent in figure 9(b). With heat- 
conduction damping, the effective circuit is that shown in figure 9(c). 

Examination of figure 8(a) shows that the maximum value of will not exceed 

0. 0075 if the maximum value of £ is made to occur at 78 percent burn time or before. 
The following conditions at 78 percent burn time were obtained by interpolation from 
tables I and n for the oxygen system: 

o> n = 102. 5 radians/sec 9 = -124° 

P g = 43.7 psi X’ = 0. 048 lb-sec/in. 5 

B u = -71 in. 5 /lb-sec r’ = 0. 083 lb-sec/in. 5 
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(a) Undamped accumulator, actual circuit. 

o — W — 

L a -j-C a 

<b) Undamped accumulator, equivalent circuit 4 
R 

° WTH1 i — VW' — | 

L 3 C 2 ^C 3 


(c) Damped accumulator. 

Figure 9. - Circuits representing accumulator 
with bellows; L 1( duct inertance; bellows 
inertance; U « + L 2 ; Ci, adiabatic com- 
pliance: C 2 , oel lows compliance; 

C a = llLd/Cj) + (1/C 2 0; Ci + C 3 , isothermal 
gas compliance; RCo = T h (thermal time con- 
stant). 


The condition that £ = £ is ^ = -1/2 0 = 62°, and the total required susceptance 

max 

B is 


B = 


X' - R* tan i p 



-12 in.^/lb-sec 


Therefore, the accumulator susceptance requirement is 

B = B - B = 59 in. Vlb-sec 
a u 

Since B & must be positive at all times, w a L a C a must be less than 1 at cutoff. The 
bellows stiffness is to be small compared to that of the gas; therefore, the net compliance 
is roughly proportional to 1/P 2 an d a t cutoff is not more than 1. 17 times its value at 
78 percent burn time. To provide a margin, o> 2 L a C a is set equal to 0.8 at cutoff. Then 
at 78 percent burn time, 


L C = — = 0. 31x10 4 sec 2 

Si Si n 

148^ x 1. 17 
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Figure 10. - Accumulator for liquid-oxygen system. 


B a = 


to C 
n a 


(o 2 h C 
n a a 


= 1 . 48 co C = 59 in. D /lb-sec 
n a 


Therefore, C = 0. 39 in. 5 /lb and L = 0. 8xl0~^ lb-sec^/in. 5 . 
a a 

A tentative design based on these parameters is shown in figure 10. Since the max- 
imum value of was selected as the design point, no further calculations are necessary 
for the undamped case. It is known that P falls quickly to negative values after 78 per- 

o 2 

cent burn time and remains negative thereafter because cu L C < 1 at cutoff. The 

n a a 

curve of f as a function of time for the modified system is sketched roughly in figure 8. 
o 

To evaluate the effect of damping by heat conduction, the same design was used 
except that the nitrogen volume was assumed to be packed with fine wire . The equivalent 
circuit is shown in figure 9(c). The accumulator admittance at 78 percent burn time 
(CO = 102. 5) was calculated by using three values of the time constant corresponding 
to T h ou n = 0. 5, 1.0, and 2.0. Values of admittance were found to be, respectively, 

15 + i86, 18 + i75, and 14 + i65 inch 5 per pound-second. Note that the susceptance is 
larger than for the undamped case because of the increase in effective compliance. Cal- 
culation of \p shows that is negative (because i J/ + 8 < -90°) at 78 percent burn 
time for all three values of R, indicating that the peak value occurs somewhat earlier 
where the envelope £ is lower. Of more importance is the effect of damping in 
reducing the peak value of £ by making X less than 1. Values of A calculated for 
the three cases (at 78 percent burn time, which is close enough to make them applicable 
at the peak) were 0.37, 0.34, and 0.39. Thus, if the design value of r h can be realized 
within a factor of two, the peak value of £ is less than 40 percent of that of the undamped 
system. With the damped accumulator the peak value of is approximately 0. 0017, 
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occurring at about 75 percent burn time. Because of the increased effective compliance, 

the damped accumulator frequency is slightly below the modal frequency at cutoff. For 

this reason, P rises above zero to 0.0031 at cutoff, as shown in figure 8(a). Reducing 
° _4 2 5 

the inertance L to 0.7x10 pound-second per inch would make P very close to 
zero at cutoff and increase the peak value only slightly to about 0. 022 at 77 percent burn 
time. At this time, for the unmodified fuel system is 0. 0014, making the maximum 
total in-phase feedback for the vehicle ^ about 0.0036 with the damped accumula- 

tor. (The fact that the peaks of P, Q for the damped system are nearly tangent to the un- 
modified curve is accidental. ) 

Additional checkout calculations would be necessary to investigate the effects of 
variations in the parameters . Modification of the fuel system also would be desirable as 
a precaution. A very small undamped accumulator with a working volume of only about 
2 cubic inches would make £ negative after about 50 percent burn time. 
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